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ABSTRACT 

The use of Type la SNe has thus far produced the most reliable measurement 
of the expansion history of the Universe, suggesting that ACDM offers the best 
explanation for the redshift-luminosity distribution observed in these events. But 
the analysis of other kinds of source, such as cosmic chronometers, gamma ray 
bursts, and high-;^ quasars, conflicts with this conclusion, indicating instead that 
the constant expansion rate implied by the i?h = ct Universe is a better £t to 
the data. The central difficulty with the use of Type la SNe as standard candles 
is that one must optimize three or four nuisance parameters characterizing su¬ 
pernova luminosities simultaneously with the parameters of an expansion model. 
Hence in comparing competing models, one must reduce the data independently 
for each. We carry out such a comparison of ACDM and the i?h = ct Universe, 
using the Supernova Legacy Survey (SNLS) sample of 252 SN events, and show 
that each model fits its individually reduced data very well. But since = ct has 
only one free parameter (the Hubble constant), it follows from a standard model 
selection technique that it is to be preferred over ACDM, the minimalist version 
of which has three (the Hubble constant, the scaled matter density and either the 
spatial curvature constant or the dark-energy equation-of-state parameter). We 
estimate by the Bayes Information Criterion that in a pairwise comparison, the 
likelihood of = cf is ~ 90%, compared with only ~ 10% for a minimalist form 
of ACDM, in which dark energy is simply a cosmological constant. Compared to 
Rh = ct, versions of the standard model with more elaborate parametrizations of 
dark energy are judged to be even less likely. 

Subject headings: cosmic background radiation - cosmological parameters - cosmology: 
observations - cosmology: theory - distance scale - supernovae: general 



1. Introduction 


Type la supernovae have a well-defined luminosity, permitting their use as standard 
candles under the assumption that the luminosities of nearby and distant sources are 
similarly related to color and light-curve shape. The use of such events has thus far produced 
the most reliable measurement of the expansion history of the Universe (Perlmutter et ah 
1998, 1999; Garnavich et ah 1998; Schmidt et ah 1998; Riess et al. 1998), leading to the 
discovery of dark energy. 

In recent years, several samples of Type la supernova events have been assembled, and 
several of these have been merged to form even larger compilations, such as the Union2.1 
sample (Kowalski et al. 2008; Suzuki et al. 2012), which currently contains 580 supernova 
detections, including events at redshift z > 1 observed with the Hubble Space Telescope 
(see, e.g., Kuznetsova et al. 2008; Dawson et al. 2009; Riess et al. 2011). However, while 
such merged samples offer several advantages, the fact that each sample has its own set 
of systematic and intrinsic uncertainties makes it difficult to fit a cosmological model to 
any merged sample. As we shall discuss, and as Kim (2011) and others have already 
pointed out, it is questionable whether the model parameters can be estimated by simply 
minimizing a since parameters characterizing error dispersion(s) must be estimated 
simultaneously. One commonly used method estimates error dispersion(s) by constraining 
the reduced y^, i.e., the y^ per degree of freedom, to equal unity. We shall find that the 
method of maximum likelihood estimation (MLE) is to be preferred; but in any method, 
the presence of multiple dispersion parameters is a complication. 

Fortunately, some of the homogeneous samples are themselves quite large, and perfectly 
suited to the type of comparative analysis we wish to carry out in this paper. In fact, 
almost half of all the Type la SNe in the Union2.1 compilation were derived from the single, 
homogeneous sample assembled during the first three years of the Supernova Legacy Survey 
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(SNLS; Guy et al. 2010). This catalog contains 252 high redshift Type la supernovae 
(0.15 < < 1.1). The multi-color light-curves of these SNe were measured using the 

MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), using 
repeated imaging of four one-square degree fields in four bands. The VLT, Gemini and Keck 
telescopes were used to conhrm the nature of the SNe and to measure their redshifts. Very 
importantly, since the same instruments and reduction techniques were employed for all 252 
events, it is appropriate to include a single intrinsic dispersion in the analysis of the Hubble 
Diagram (HD) constructed from this sample. The study of catalogs such as this has led 
to a general consensus that AGDM offers the best explanation for the redshift-luminosity 
relationship, and observational work is now focused primarily on rehning the fits to improve 
the precision with which the model parameters are determined. This is one of the principal 
motivations for attempting to merge samples to create catalogs with broader redshift 
coverage and better statistics. 

But as successful as this program has been, several drawbacks associated with the 
use of Type la supernovae have made it necessary to seek alternative methods of probing 
the cosmic spacetime. It is quite difficult to use supernova measurements in unbiased, 
comparative studies of competing expansion histories, since at least three or four ‘nuisance’ 
parameters characterizing the standard candle must be optimized simultaneously with each 
model’s free parameters, rendering the data compliant to the underlying theory (Kowalski 
et al. 2008; Suznki et al. 2012; Melia 2012a). Several notable attempts have been made to 
mitigate the impact of this model dependence, e.g., through the use of kinematic variables 
and geometric probes that avoid parametrizing the hts in terms of pre-assumed model 
components (Shafieloo et al. 2012). In the end, however, even these models reqnire the 
availability of measnrements based on standard candles. Unfortunately, Type la SNe may 
be nsed for this pnrpose only as long as the nnisance parameters characterizing their 
lightcurves are known. The application of such methods to the Union2.1 Type la SN sample 
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uses nuisance parameters optimized for the concordance model, so the current results are 
not completely free of any biases. 

The expansion of the Universe is now being studied by several other methods, including 
the use of cosmic chronometers (Jimenez & Loeb 2002; Simon et ah 2005; Stern et ah 2010; 
Moresco et ah 2012; Melia & Maier 2013), gamma ray bursts (GRBs; Norris et ah 2000; 
Amati et ah 2002; Schaefer 2003; Wei & Gao 2003; Yonetoku et ah 2004; Ghirlanda et ah 
2004; Liang & Zhang 2005; Liang et al. 2008; Wang et al. 2011; Wei et ah 2013), and high-;^ 
quasars (Kauffmann & Haehnelt 2000; Wyithe & Loeb 2003; Hopkins et al. 2005; Groton 
et al. 2006; Fan 2006; Willott et al. 2003; Jiang et al. 2007; Kurk et al. 2007; Tanaka & 
Haiman 2009; Lippai et al. 2009; Hirschmann et al. 2010; Melia 2013a). In contrast to the 
perception based on Type la SNe that AGDM can best account for the observed expansion 
of the Universe, the conclusion from these other studies is that the cosmic dynamics is 
better described by a cosmology we refer to as the = ct Universe (Melia 2007; Melia & 
Abdelqader 2009; Melia & Shevchuk 2012). 

For example, a comparative analysis was recently carried out of the AGDM and 
= ct cosmologies using the GRB Hubble diagram (Wei et al. 2013). This study found 
that once the various parameters are estimated for each model individually, the i?h = ct 
cosmology provides a better £t to the data. And although about 20% of the GRB events 
he at least 2a away from the best-fit curves (in both models), snggesting either that some 
contamination by non-standard GRB Inminosities is nnavoidable, or that the errors and 
intrinsic scatter associated with the data are being underestimated, various model selection 
techniqnes applied to the GRB data show that the likelihood of = ct being the correct 
model rather than AGDM is ~ 90% versus ~ 10%. 

Moreover, although the redshift-distance relationship is essentially the same in AGDM 
and the i?h = ct Universe (even out to 2 : > 6 — 7), the redshift-age relationship is not. 
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This motivated a recent examination of whether or not the observed growth rate of high-;^ 
quasars could be used to test the various models. Quasars at z > 6 are now known to be 
accreting at, or near, their Eddington limit (see, e.g., Willott et ah 2010a,b; De Rosa et ah 
2011), which presents a problem for ACDM because this makes it difficult to understand 
how ~ 10® Mq supermassive black holes could have appeared only 700-900 Myr after the 
big bang. Instead, in i?h = ct, their emergence at redshift ~ 6 corresponds to a cosmic age 
of > 1.6 Gyr, which was enough time for them to begin growing from ~ 5 — 20 Mq seeds 
(presumably the remnants of Pop II and III supernovae) at < 15 (i.e., after the onset of 
re-ionization) and still reach a billion solar masses hy z ^ 6 via standard. Eddington-limited 
accretion (Melia 2013a). 

In light of this apparent conflict between the implications of the Type la supernova 
work and the results of other studies, we have begun to look more closely at the possibility 
of directly comparing how ACDM and i?h = ct account for the supernova measurements 
themselves. This is not an easy task, principally because of the enormous amount of work 
that goes into hrst establishing the SN magnitudes, and then carefully htting the data using 
the comprehensive set of parameters available to ACDM (see, e.g., Suzuki et ah 2012). In a 
previous paper (Melia 2012a), the redshift-distance relationship in ACDM was compared 
with that predicted by = ct, and each with the Union2.1 sample. It was shown that the 
two theories produce virtually indistinguishable prohles, though the £t with = ct had not 
yet been fully optimized. That is, despite the fact that this previous analysis simply used 
Rh = ct to £t the data optimized with ACDM, the results were quite promising, suggesting 
that a full optimization procedure—followed separately for = ct and ACDM—ought to 
be carried out. The principal goal of this paper is to complete this study. 

A direct comparison of ACDM with = ct using Type la SNe is also motivated by 
recent theoretical work suggesting that the Friedmann-Robertson-Walker (FRW) metric 
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is more specialized than was previously thought (Melia 2013b). A close examination of 
the physics behind the symmetries incorporated into this well-known and often employed 
solution to Einstein’s equations showed that the FRW metric applies only to a fluid with 
zero active mass, i.e., a fluid in which p -|- 3p = 0, where p is the total pressure and p the 
total energy density. This is consistent with the equation-of-state used in = ct, but not 
in ACDM, except that when one averages the pressure, (p), over the age of the Universe, 
one does get ~ — (p)/3 for the estimated ACDM parameters. ACDM therefore appears to 
be an empirical approximation to i?h = cf, asymptotically approaching the requirements of 
the zero active mass condition consistent with the FRW metric. 

In § 2 of this paper, we briefly summarize the contents of the SNLS sample, and explain 
the model parameters (including those associated with the data) that are to be estimated. 
We present the fits to the supernova data in § 3, and a direct comparison between ACDM 
and the Rh = cf Universe is made in § 4. We end with a discussion and conclusions in § 5. 


2. The SNLS Supernova Sample 

The Supernova Legacy Survey (SNLS) sample contains 252 high-redshift 
(0.15 < 2 : < 1.1) Type la SNe discovered during the hrst three years of operation 
(Guy et ah 2010). One of the most important features of this catalog is that it constitutes 
a single, homogeneous sample. It also covers the very important redshift range where 
the standard model suggests the Universe underwent a transition from decelerated to 
accelerated expansion. Guy et al. (2010) used two light-curve htters (SALT2 and SiFTO) 
to determine the peak magnitudes, light-curve shapes, and colors of the Type la SNe. 

For SiFTO, a distance modulus is dehned for each SN as the linear combination 

= rriB + a ■ {s - 1) - I3 ■ C - Mb, 


( 1 ) 



where is the peak rest-frame 5-band magnitude, s is the stretch (a measure of 
light-curve shape), C is the color (peak rest-frame B — V), and Mb is the absolute 
magnitude of a Type la supernovalU When corrected for shape and color. Type la SN 
luminosities have a dispersion of only ~ 15%. The coefficients a,/9 are thus the parameters 
of a luminosity model, though it is not a full statistical model, since it lacks an explicit error 
dispersion parameter. In the present context, a, (3 and Mb are ‘nuisance’ parameters, as 
they cannot be estimated independently of an assumed cosmology. They must be optimized 
simultaneously with the cosmological parameters, as will be explained in § 3. 


The theoretical distance modulus /xth is calculated for each SN from its measured 
redshift z by 

hth(^) = 51og(T>L(2;)/10pc), (2) 

where Dl{z) is the mo del-dependent luminosity distance. A determination of Dl requires 
the assumption of a particular expansion scenario. Both ACDM and 5h = ct are FRW 
cosmologies, but the former assumes specific constituents in the density, written as p = 

Pr + Pm + Pa, where Pr,Pm,PA are, respectively, the energy densities of radiation, matter 
(luminous and dark), and the cosmological constant A. These are often expressed in terms of 
today’s critical density, pc = Sc^Hq/StiG, where Hq is the Hubble constant, by Qm = Pm/Pc, 
hlr = Pr/Pc, and IIa = Pa/P c- In a flat universe with zero spatial curvature, the total scaled 
energy density is H = Vtm + fir + Ha = 1- Since 1, -|- Ha = 1. In = ct, on 

the other hand, whatever constituents are present in p, the principal constraint is the total 
equation-of-state p = —p/3, which as we mentioned in § 1, is in fact required by the use of 


^The SNLS supernova compilation of 252 SNe is currently available in the University of 
Toronto’s Research Repository, at https://tspace.library.utoronto.ca/handle/1807/24512. It 
includes the following information for each SN: ttib (with corresponding standard error (Tm^), 
s (with error cXs), C (with error ac), and the covariances between mB,s,C. 





9 


the FRW metric. 


When dark energy is included with an unknown equation-of-state, p\ = w\px, the 
general form of the luminosity distance in ACDM is given by 


^ACDM(;^) 


c {1 + z) 

^0 a/I fife 


smn< 



X 


dz 

V^l + + ff/l + x)2 + + ;2)3(1 +-a) 

where c is the speed of light. In this, hlfc = 1 — flm — IIa represents the spatial curvature of 
the Universe—appearing as a term proportional to the spatial curvature constant k in the 
Friedmann equation. In addition, sinn is sinh when hlfc > 0 and sin when Qk < 0. For a flat 
Universe {Qk = 0), the right side becomes (1 + z)c/Hq times the indehnite integral. 




In the i?h = ct Universe, the luminosity distance is given by the much simpler 
expression 

= A(1 + ^) HI+^). (4) 

-no 

The factor c/Hq is in fact the gravitational horizon i?h(^o) (which itself is coincident with 
the Hubble radius) at the present time, so we may also write the luminosity distance as 

D^^='^\z) = Rh(to)(l + z) ln(l + z) . (5) 


A more extensive description of the observational differences between ACDM and i?h = ct 
is provided in Melia (2007,2012a), Melia & Shevchuk (2012), Melia & Maier (2013), and 
Wei et al. (2013). For a pedagogical treatment, see also Melia (2012b). 


3. Theoretical Fits 

In this paper, we shall use two methods for calculating point and interval estimates 
of cosmological model parameters: one a method commonly used in past analyses (e.g.. 
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Kowalski et al. 2008; Guy et al. 2010; Suzuki et al. 2012), and one more recently proposed, 
which is based on maximizing the likelihood function (e.g., Kim 2011; Melia & Maier 2013; 
Wei et al. 2013). As we shall see, the latter allows us to estimate all parameters, including 
the unknown intrinsic dispersion. But we shall hnd that under some circumstances, 
the manner in which this latter method works offers some justihcation for the former, 
which instead, minimizes subject to the condition that the reduced equal unity. In 
either method, however, nuisance parameters characterizing the SN luminosities, and free 
parameters of the cosmological model, must be htted simultaneously. 

Initial Remark on Free Parameters. Before we begin the actual htting of the various 
models, we take a moment to summarize the free parameters available in each. As alluded 
to earlier, it is now understood that the symmetries incorporated into the FRW metric 
require a zero active mass condition (i.e., p + 3p = 0) in the cosmic fluid (Melia 2013b). 
One can understand this even without a formal proof, by using the following reasoning. 

In spherically symmetric spacetimes, which include the special FRW case, a proper 
mass emerges from the introduction of the metric into Einstein’s held equations. We call this 
the Misner-Sharp mass (Misner & Sharp 1964), which is dehned as M(i?h) = (47r/3)i?j^p/c^, 
in terms of the gravitational radius i?h = ‘2GM{R\^/c^. This proper mass is a consequence 
of Birkhoff’s theorem (Birkhoff 1923), which states that none of the mass-energy beyond i?h 
contributes to the spacetime curvature within a shell at this radius. The gravitational 
radius must therefore be dehned as written here, although the cosmos itself may be inhnite. 
But since M is a proper mass dehned in terms of the proper density p and proper volume, 
Rh must therefore be a proper radius as well. Thus, by Weyl’s postulate (Weyl 1923), this 
gravitational radius must have the form Rh = where a{t) is the universal expansion 
factor and rh is a constant comoving distance. One can therefore easily show from the hrst 
Friedmann equation (see, e.g., Melia & Shevchuk 2012) that a? + kc^ = r^'^, where k is 



-li¬ 
the spatial curvature constant. In other words, the expansion rate a must be constant, 
which then leads, via the second Friedmann (or ‘acceleration’) equation, to the condition 
p = —p/3. The equation-of-state w = —1/3 is therefore required in this cosmology, and 
cannot be adjusted to optimize hts to the data. 

Further, since H = (d/a) = 1/t, one can easily see that the Hubble radius 
.Rnubbie = c/H must be given as i?Hubbie = ct in this cosmology. But in general relativity, we 
recognize such a surface moving at proper speed c relative to an observer in the co-moving 
frame as an event horizon. So if the Misner-Sharp mass gives rise to the gravitational 
radius which is identihed with i?Hubbie, then from the Friedmann equation we must also 
have k = 0. The = ct Universe is therefore flat and has a unique, hxed equation-of-state 
w = —1/3. It has only one free parameter: the Hubble constant Hq. 

In ACDM there is some flexibility in choosing our free parameters for the analysis of 
Type la SNe, depending on how we characterize dark energy. At the very minimum, we may 
adjust Hm, k (or, equivalently, Hde), Wde and Hq. (Note that additional free parameters, such 
as the baryon density, may be relevant to hts involving the CMB, but are not necessary for 
constructing the SN Hubble diagram.) And though inhation is not yet a fully self-consistent 
model, we may include it as a foundation for ACDM, in which case we could reasonably 
argue that k = 0. However, even inhation does not specify the equation-of-state for dark 
energy. And using priors from other kinds of observations is not really appropriate for this 
work, because their values change from observation to observation—and from instrument to 
instrument. For example, note that the values of Wde and Dm from WMAP (Bennett et ah 
2003) are quite diherent from those inferred by Planck (Ade et al. 2014). The most basic 
ACDM model we may use for supernova work thus has three free parameters. Below, we 
shall choose the parameters Dm, tCde and Hq to represent the standard model, since these 
are the three that produce the most favorable hts when using ACDM, based on what is 
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currently known. 


Method I. The commonly employed fitting method, as discussed by Guy et ah (2010), 
determines the best-fit cosmology by (iteratively) minimizing the function 




h'th(^i)) 
^lc,i + ^int 


( 6 ) 


over a joint parameter space, thereby optimizing both the nuisance parameters 
and the cosmological parameters. Here i ranges over the serial numbers of the SNe in the 
SNLS catalog, and depends on a, (3, Mb according to Equation 1; and of course, 
depends on the cosmological parameters. The SN-speci£c dispersion ^ is defined by 


a 


lc,2 


2 I 2 2 

= 


+ 


2^2 

^C,i 


+ a 


ms sC,i ? 


(7) 


where crmgx o's,i, and ae,i are the standard errors of the peak magnitude and light-curve 
parameters of the SN. The term Cmgsc,! comes from the covariances among mB,s,C, 
and likewise depends quadratically on a, [3. The quantity Uint is an intrinsic (i.e., 
SN-independent) dispersion, the value of which is set by requiring that the reduced equal 
unity. This additional dispersion takes into account, e.g., the residual (~ 15%) variability 
in the Type la SN luminosity, not captured by the correction coefficients a, [3. 

The following is how Method I is applied, yielding fitted values for the cosmological 
parameters (with the exception of the Hubble constant Hq —more on this below), and 
model-specific fitted values for the nuisance parameters a, (3, Mb- First, Uint is initialized to 
zero. Second, the function is minimized over the cosmological and nuisance parameters, 
and their values are updated. Third, with these updated values fixed, the value of Uint for 
which the reduced i.e., y^/dof = x^/(n — k), equals unity is determined, and a-mt is 
updated. (Here, n [= 234, see below] is the number of SNe and k is the total number 
of parameters, including Uint-) Steps 2 and 3 are repeated until the parameters converge 
within tolerance to stable values. 
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To find the best-fit coefficients a, /5, Mb and the cosmological parameters that dehne 
the fitted model, we use Markov-chain Monte Carlo (MCMC) techniques in our calculations. 
Our MCMC approach generates a chain of sample points distributed in the parameter 
space according to the posterior probability, using the Metropolis-Hastings algorithm 
with uniform prior probability distributions, such that 0.5 < a < 2.0, 2.0 < /3 < 4.0, 

— 19.5 < Mb < —18.5, 0.0 < (Tint < l-oj^ 0.0 < Om < 1-0, and/or 0.0 < Oa < 1.8 (in 
the non-flat case). In the parameter space formed by the constrained nuisance parameters 
and cosmological parameters, a random set of initial values of the model parameters is 
chosen to calculate the {Method I), or the likelihood function {Method II). Whether the 
set of parameters can be accepted as an effective Markov chain or not is determined by 
the Metropolis-Hastings algorithm. The accepted set not only forms a Markov chain, but 
also provides a starting point for the next process. We then repeat this process until the 
established convergence accuracy can be satisfied. For each Markov chain, we generate 10^ 
samples consistent with the likelihood function. Then we derive the probability distributions 
of the coefficients a, (I, Mb, and the cosmological parameters from a statistical analysis of 
the samples. 

To be consistent with previous work (Guy et ah 2010), we discard from the SNLS 
catalog all SNe with a peak rest-frame {B — V) > 0.2. Such red SNe are found only at 
2 ; < 0.6 in the SNLS because they are fainter than the average, and hence are undetected (or 
unidentihed spectroscopically) at higher redshifts. Discarding them minimizes any potential 
biasing of the distance modulus by an inadequate color correction. Indeed, the correction 
coefficient (/3) we estimate from the bulk of the SNe may not apply to those red SNe that 
are more likely to be extinguished by dust in their host galaxy than bluer ones. This cut, 

^In Method II, (Tint is treated on the same level as the other parameters—more on this 


below. 
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applied to both SALT2 and SiFTO samples, discards 11 SNe. There are also 3 SNe, the 
peak magnitudes of which could not be obtained with SiFTO due to a lack of observations 
in the qm and zm bands. Finally, again following Guy et ah (2010), we remove SNe that 
are 3(T outliers for either of the best-£t ACDM and i?h = ct models. This removes 4 more 
SNe from both the ACDM and i?h = ct samples. Three (o4D3dd, 04Dlak, and 04D4gz) are 
in common; 03D4cx is removed for ACDM, while 05D2ei is removed for R\^ = ct. In total, 
that leaves n = 234 SNLS supernovae for the analysis. 

Although the number of free parameters in the dark-energy model can be as large as 
eight, depending on how one handles the dark energy, in this paper we take a minimalist 
approach and use only three of the most essential ones; these are taken from among a 
set that includes the Hubble constant the matter energy density Dm, the dark energy 
density Da, and the dark energy equation-of-state parameter wj^. In fitting, the chosen 
value of Hq is not independent of Mb- That is, one can vary either of Hq, Mb, but not both. 
Therefore, if we take Mb to be a member of the set of nuisance parameters characterizing 
the SN luminosities, the dark-energy model will have, at most, only two parameters: Dm 
and Da. For consistency with Guy et al. (2010), we use Hq = 70 km s“^ Mpc“^ (recalling 
that this value does not actually affect any of the fits). 

To ensure that our implementation of Method I is reliable and consistent with the 
results of Guy et al. (2010), we first examine a simple test case: the flat Universe, in 
which Da = 1 — Dm, with equation-of-state parameter = — 1- Dm is then the only free 
cosmological parameter. Guy et al. (2010) found that the assumption of flatness resulted 
in a value for Dm that we now know is inconsistent with that reported by Planck2013 (Ade 
et al. 2014). Nonetheless, it is useful to begin with this test case because it allows us to 
compare our results directly with theirs. Applying the above y^-minimization procedure, 
we find that the fitted parameters in this case are a = 1.470 ± 0.121, (3 = 3.207 ± 0.157, 
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Fig. 1.— The (normalized) likelihood distributions and 2-D joint distributions with la, 2cr 
contours, for the coefficients a, (5, Mb in Equation 1 and the matter density parameter hlnn 
obtained from a constrained minimization of {Method I), with the assumption of flatness 
(IIa = 1 — him)- The la confidence intervals are well approximated by Gaussians. 
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Fig. 2.— The (normalized) likelihood distributions and 2-D joint distributions with la, 2a 
contours, for the coefficients a, f3, Mb in Equation 1 and the cosmological parameters him, 
obtained from a constrained minimization of {Method I), without the assumption of 
flatness. 
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Mb = —19.197 ± 0.031, and f2m = 0.235 ± 0.035, with dint = 0.090. In Figure 1 we show 
the (normalized) likelihood distribution for each parameter (a, (3, Mb] flm), according to the 
factor exp(—x^/2), and la,2a contours for the joint distribution of each pair of parameters. 
For each parameter, the likelihood distribution is well approximated by a Gaussian, and the 
stated conhdence interval is a 68% (i.e., ila) interval for this Gaussian. Using the same 
method and conditions, Guy et ah (2010) obtained a = 1.487 ± 0.097, j3 = 3.212 ± 0.128, 
Mb = —19.210 ± 0.030, and Gm = 0.215 ± 0.033, with Uint = 0.087. These results are quite 
consistent, again attesting to the reliability of our calculation. 

To allow for a more realistic £t to the SNLS supernova data, we removed the flatness 
assumption, and allowed both Gm and Ga to be free parameters. By minimizing we now 
obtain a = 1.469±0.123, (3 = 3.209±0.159, Mb = -19.187±0.068, and = 0.217±0.150, 
Ga = 0.718 ± 0.329, with dint = 0.090. The (normalized) likelihood distribution for each 
parameter (a, /9, Gm, Ga), and a contour plot for each two-parameter combination, are 
shown in Figure 2. Dropping the flatness condition, i.e., fitting both and Da, makes the 
best-fit AGDM model marginally consistent with the Planck2013 results, because the la 
standard errors are now considerably larger than in the previous case. 

It should be noted that the approach of Method I does not actually maximize 
a likelihood. Rather, it calculates a value for Uint (with no accompanying uncertainty) 
by requiring that x\oi equal unity. Driving to unity, irrespective of how well the 
cosmological model fits the data, makes it impossible to perform a fair comparison of 
competing models by using Method I, especially when they have different numbers of 
parameters. A statistically valid analysis must estimate error dispersion parameter(s), along 
with all other parameters, by maximizing a joint likelihood function. 

Method II. As a superior alternative to Model I for estimating cosmological parameters, 
and (simultaneously) the model-specific optimized nuisance parameters a, (3, Mb, we 
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employ a method described in D’Agostini (2005) and Kim (2011), which is based on the 
maximization of a joint likelihood function. The joint likelihood function for all these 
parameters and the intrinsic dispersion (Tint, based on a flat Bayesian prior, is 


L=n 


^27r(a2 . + 


X exp 


(h'Bji h'th(^i)) 
+ <^int) 


oc exp (-xV2) • 


( 8 ) 


As in Method /, each distance modulus depends on a, (3, Mb, and the theoretical 
distance modulus iJiih{,Zi) depends on the cosmological parameters. Maximum likelihood 
estimation (MLE) is, of course, a standard statistical procedure, and appeared in Method I 
in a limited way, as a minimization of The new feature in Method II is that (Tint is 
treated on the same level as the other parameters: the likelihood L is maximized over all 
parameters, now including (Tint- This ‘full MLE’ provides a statistically founded method 
for estimating (Tint (Kim 2011). It also treats on the same level the cosmological parameter 
uncertainties and the potential uncertainty in (Tint, which can affect each other. 


The optimized ACDM parameters obtained by Method II are: a = 1.275 ± 0.120, 

(3 = 2.637T0.155, Mb = -19.165±0.081, = 0.365T0.137, and IIa = 0.846T0.353, with 

o'int = 0.103 ± 0.010. For each parameter, the likelihood distribution is well approximated 
by a Gaussian, and the stated conhdence interval is a 68% (i.e., ±1(t) interval for this 
Gaussian. As anticipated by Kim (2011), these best-fit values are not all consistent with 
those of Method /, though the estimated value of Gm and its error are entirely consistent 
with the Planck2013 results. In an analysis using mock samples, Kim concluded that 
differences such as these arise because the commonly used method based on the constrained 
minimization of (i.e.. Method I) does not include in its error propagation the covariance 
of (Tint with the other parameters. 


The Method II values of the coefficients a, (3, Mb may now be used to calculate the 
distance modulus /is of Equation 1, for each SN. From the distance moduli, we construct 
the Hubble diagram shown in the upper panel of Figure 3. We also show the corresponding 
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SNLS sample of 234 Type la SNe. As is now well known, the theoretical ht is excellent. 
The maximum value of the joint likelihood function is given by —21nL = —238.40, which 
we shall need when comparing models using the Bayes Information Criterion (see below). 
For completeness. Figure 3 also shows the Hubble diagram residuals corresponding to 
the best-£t ACDM model. The likelihood distribution obtained by Method II for each 
parameter (a,/9, Hm, Ha; o'int), and a contour plot of the joint distribution for each 
two-parameter combination, are shown in the lower panel. 


4. A Direct Comparison Between ACDM and the = ct Universe 


In the i?h = ct Universe, there is only one free parameter: the Hubble constant Hq. But 
since we cannot vary Mb and Hq separately, there are no free parameters left to adjust the 
theoretical curve, once we optimize Mb as one of the nuisance parameters characterizing 
the supernova luminosities. 

Using the MLE approach {Method II), we find that for the i?h = ct Universe, 
the optimized nuisance parameters are a = 1.175 ± 0.115, fl = 2.608 ± 0.149, 

Mb = —18.959 ± 0.011, with dint = 0.106 ± 0.010. All likelihood distributions are well 
approximated by Gaussians, and the given confidence intervals are 68% (i.e., ±la) intervals 
for the Gaussians. As with ACDM, we plot the (normalized) likelihood distribution for 


Table 1. Best Fits for Different Cosmological Models 


Model 

a 

fs 

Mb 

c^int 

Qm 

Oa 

-2InL 

BIC 

Method I 









ACDM (fc = 0) 

1.470 ± 0.121 

3.207 ± 0.157 

-19.197 ± 0.031 

0.090 

0.235 ± 0.035 

1 - Qni 

— 

— 

ACDM 

1.469 ± 0.123 

3.209 ± 0.159 

-19.187 ± 0.068 

0.090 

0.217 ± 0.150 

0.718 ± 0.329 

— 

— 

Method II 









ACDM 

1.275 ± 0.120 

2.637 ± 0.155 

-19.165 ± 0.081 

0.103 ± 0.010 

0.365 ± 0.137 

0.846 ± 0.353 

-238.40 

-205.67 


1.175 ± 0.115 

2.608 ± 0.149 

-18.959 ± 0.011 

0.106 ± 0.010 

— 

— 

-231.85 

-210.03 







0.2 0.4 0.6 0.8 1.0 



Fig. 3.— Upper panel: Hubble diagram and Hubble diagram residuals for the SNLS sample 
of 234 Type la SNe in ACDM. Solid curve: best-fit model with = 0.365, Ha = 0.846 and 
Pa = —Pa- Lower panel: (normalized) likelihood distributions and 2-D joint distributions 
with la, 2a contours, for the coefficients a, (3, Mb, the intrinsic dispersion cxint, and Hin,HA. 
The htting method employed is MLE {Method II). 
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45 -1-1-1-1-1-1-1-1-r 



-0.8 -1-'-1-■-1-■-1-■-L_ 

0.2 0.4 0.6 0.8 1.0 



Fig. 4.— Same as Figure 3, but now for the = ct Universe. Upper panel; Hubble diagram 
and Hubble diagram residuals for the data optimized for = ct, and the corresponding 
theoretical curve. Lower panel; The (normalized) likelihood distributions and 2-D joint 
distributions with Icr, 2a contours, for all fitted parameters. The fitting method employed is 
MLE {Method II). 
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each parameter (a,/9, (Tint), and 2-D plots for two-parameter combinations (Fignre 4). 
The best-£t valnes are qnite similar to those for ACDM, but are not exactly the same, 
reaffirming the importance of reducing the data separately for each model being tested. The 
= ct distance modulus is compared to the SNLS sample in the upper panel of Figure 4. 
For completeness, we also show the Hubble diagram residuals corresponding to the = ct 
Universe at the bottom of this panel. The maximum value of the joint likelihood function 
for the optimized i?h = ct fit corresponds to —21nL = —231.85. All the fits performed in 
this paper are summarized in Table 1, for ease of comparison. 

An inspection of the Hubble diagrams in Figures 3 and 4 reveals that the distance 
moduli are slightly different when the nuisance parameters are optimized using different 
models, but both ACDM and = ct fit their respective data sets very well. One certainly 
gets this impression from a side-by-side comparison of the Hubble-diagram residuals for 
the SNLS sample in ACDM and i?h = ct, shown in Figure 5. However, because these 
models formulate their observables (such as the luminosity distance in Equations 3 and 5) 
differently, and because they do not have the same number of free parameters, a decision 
between the models must be based on a formal model selection technique, and in this 
regard, the results of our analysis favor i?h = ct over ACDM, as we shall now demonstrate 
quantitatively. 

A companion paper (Melia & Maier 2013) discussed at length how one may use 
state-of-the art model selection tools to choose the model to be preferred in accounting 
for the data. We shall not reproduce that discussion here, but we do point out that to 
assess competing models in cosmology, a strong case for using the Bayes Information 
Criterion (BIC) has been made (see, e.g.. Fiddle 2004, 2007; Fiddle et al. 2006). The BIC 
is applicable when data points are independent and identically distributed, which is a 
reasonable assumption for supernova redshift-luminosity data. The method has now been 
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used to compare several popular models against ACDM (see, e.g., Shi et al. 2012). 



Fig. 5.— A side-by-side comparison of the Hubble-diagram residuals for the SNLS sample 
in the standard (ACDM) model and the i?h = ct Universe. 

Despite its name, the BIC is based not on information theory, but rather on an 
asymptotic (n —)■ cxd, where n is the number of data points) approximation to the outcome 
of a Bayesian inference procedure for deciding between competing models (Schwarz 1978). 
A comparison between the BIC values of two or more htted models provides their relative 
ranks, and also a numerical measure of conhdence (a likelihood or posterior probability) 
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that each model is the best. Unlike some model selection techniques, the BIG can be applied 
to ‘non-nested’ models, such as those we have here. The BIG for a htted regression model, 
linear or nonlinear, captures the dominant (in the n —)■ cxd limit) behavior of the associated 
Bayes factor, which assesses the strength of the evidence in its favor (Kass & Raftery 1995, 
§4.1). In this limit, it is only subdominant terms that are affected by such considerations 
as the choice of a Bayesian prior (Kuha 2004), or the extent of the model nonlinearity if any 
(Haughton 1988). There is accordingly a simple formula for the BIG, namely 

BIG =-21nL+(Inn)A; , (9) 

where k is the number of htted parameters. The logarithmic penalty term in the BIG 
strongly suppresses overhtting if n is large (the situation we have here, with n = 234, which 
is deep in the asymptotic regime). 


A quantitative ranking of htted models 1 and 2 is computed as follows. If BIG^ 

comes from model a, the unnormalized conhdence in this model is the ‘Bayes weight’ 

exp(—BIG q/ 2). That is, in the light of the data, model a has likelihood 

p/ X _ exp(-BIG^/2) _ 

exp(—BIGi/2) + exp(—BIG 2 / 2 ) 

of being the correct choice. The strength of the evidence for model 1 and against model 2 is 
quantihed by ABIG = BIG 2 — BIGi, and the following qualitative interpretation of ABIG 
is standard. If ABIG is less than 2, the evidence is “not worth more than a bare mention” 
(Kass & Raftery 1995). If it is in the range 2... 6, the evidence is positive; if it is in the 
range 6 ... 10, the evidence is strong; and if it is greater than 10, the evidence is very strong. 


With n = 234 data points and k = A parameters, the BIG for the optimized i?h = ct 
Universe is BIGi = —210.03. For the optimized AGDM, with k = 6, the corresponding 
value is BIG 2 = —205.67. Our analysis therefore shows that the evidence supplied by the 
SNLS sample for the i?h = ct Universe over AGDM is positive; and quantitatively, i?h = ct 
is favored over AGDM with a likelihood of ~ 90% versus only ~ 10%. 
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5. Discussion and Conclusions 

Our comparative analysis of ACDM and the i?h = ct Universe using the SNLS sample 
has shown that—contrary to earlier claims—the Type la SNe do not point to an expansion 
history of the Universe in conflict with that implied by other kinds of sonrce, such as the 
cosmic chronometers (Jimenez & Loeb 2002; Simon et ah 2005; Stern et al. 2010; Moresco 
et ah 2012; Melia & Maier 2013), GRBs (Norris et ah 2000; Amati et al. 2002; Schaefer 
2003; Wei & Gao 2003; Yonetoku et al. 2004; Ghirlanda et al. 2004; Liang & Zhang 2005; 
Liang et al. 2008; Wang et al. 2011; Wei et al. 2013), and high- 2 : qnasars (Kanffmann & 
Haehnelt 2000; Wyithe & Loeb 2003; Hopkins et al. 2005; Groton et al. 2006; Fan 2006; 
Willott et al. 2003; Jiang et al. 2007; Knrk et al. 2007; Tanaka & Haiman 2009; Lippai et al. 
2009; Hirschmann et al. 2010; Melia 2013a). A difficulty with the use of supernova data has 
always been their dependence on the nnderlying cosmology. There is no question that to 
compare different models properly, one must optimize the nuisance parameters describing 
the snpernova luminosities separately for each expansion scenario. It is not appropriate to 
use data optimized for AGDM to test other models. 

This has been the primary focus of our analysis in this paper. We have confirmed 
the argnment made by Kim (2011), in particular, that one should optimize parameters 
by carrying out a maximum likelihood estimation in any situation where the parameters 
include an nnknown intrinsic dispersion. The commonly used method, which estimates the 
dispersion by requiring the reduced to equal unity, does not take into account all possible 
covariances among the parameters. In this regard, our best-fit models for AGDM do not 
agree exactly with those of Guy et al. (2010), who optimized the model parameters using 
only Method I. Indeed, while their best-fit model is characterized by parameters noticeably 
different from those of Planck2013, we have demonstrated that the use of Method II actually 
results in an optimized AGDM model consistent with that based on the analysis of the 
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CMB fluctuations. Simply based on a consideration of the standard model, therefore, our 
results support the proposal made by Kim (2011) that Method II should be preferred over 
Method I in the analysis of Type la SN samples that include unknown intrinsic dispersions. 

More importantly, we have found that, when the parameter optimization is handled 
via the joint likelihood function, both ACDM and i?h = ct ht their individually optimized 
data very well. However, the i?h = ct Universe has only one free parameter—the 
Hubble constant Hq —which enhances the signihcance of the £t over that using ACDM. 
Indeed, standard model selection techniques penalize models with a large number of free 
parameters. We have found that the BIC favors i?h = ct over ACDM, by a likelihood 
of ~ 90% versus ~ 10%. The difference in likelihoods would be even greater for other 
variations of ACDM that include a larger number of free parameters, e.g., to characterize 
the dark-energy equation-of-state, should wa not be constant. 

This result would be quite significant on its own. But when we consider it in concert 
with the analysis of all the other data sets that have been analyzed thus far, it is reasonable 
to conclude that i?h = cf is a more accurate representation of the Universe than is ACDM. 
This has far-reaching consequences that will be addressed at greater length elsewhere. 

Obviously, these results call into question the conclusion that the Universe is currently 
undergoing a period of acceleration, following an earlier period of deceleration. The fact 
that the fits to the data using ACDM often come very close to those of i?h = ct (see 
Figure 5) lends weight to our suspicion that the standard model functions as an empirical 
approximation to the latter, since it has more free parameters and lacks that essential 
ingredient in = ct. the equation-of-state p = —p/3. In attempting to identify the reasons 
why ACDM produces phases of acceleration and deceleration, one is reminded of trying to 
fit a straight line with a low-order polynomial—it is always possible to make the ends meet, 
but there will be inevitable wiggles in between. It appears that the early deceleration and 
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current acceleration indicated by ACDM are two of these wiggles, whereas = ct hts the 
straight line perfectly. 

In this work, we have also demonstrated that the class of Type la SNe continues to be 
critical to our understanding of how the Universe evolves, the results reported here affirming 
the expectation from theory and general relativity that only a perfect fluid with zero active 
mass (i.e., p + 3p = 0) can be consistent with the use of an FRW metric (Melia 2013b). 
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